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Abstract. We consider a collisionless ensemble of classical particles coupled 
to a Klein-Gordon field. For the resulting nonlinear system of partial differ- 
ential equations, the relativistic Vlasov-Klein-Gordon system, we prove local- 
in-time existence of classical solutions and a continuation criterion which says 
that a solution can blow up only if the particle momenta become large. We 
also show that classical solutions are global in time in the one-dimensional 
case. 



1. Introduction 

In kinetic theory one often considers collisionless ensembles of classical parti- 
cles which interact only by fields which they create collectively. This situation is 
commonly referred to as the mean field limit of a many particle system. Such sys- 
tems have been studied extensively. In the case of non-relativistic, gravitational or 
electrostatic fields the corresponding system of partial differential equations is the 
Vlasov-Poisson system, in the case of relativistic electrodynamics it is the Vlasov- 
Maxwell system and in the case of general relativistic gravity the Vlasov-Einstein 
system. 

On the other hand the coupling of a single particle to a classical or quantum 
field has been studied. In case of the Maxwell field this is a classical problem pQ, 
but the actual dynamics and asymptotics of such systems is still an active area 
of research JH] \7\ El |5| ■ In [S] the case of a single classical particle coupled to a 
quantum mechanical Klein-Gordon field was investigated. 

In the present paper we consider a collisionless ensemble of particles moving at 
relativistic speeds, coupled to a Klein-Gordon field. This is a natural generalization 
of the one-particle situation just described. Let / = f(t,x,v) > denote the density 
of the particles in phase space, p = p(t,x) their density in space, and u — u(t,x) 
a scalar Klein-Gordon field; tgR, i€R 3 , and vSK 3 denote time, position, and 
momentum, respectively. The system then reads as follows: 

d t f+v-d x f-d x u-d v f = 0, (1.1) 



d 2 t u~ Au + u = -p, (1.2) 
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p(t,x) = j f(t,x,v)dv. (1.3) 

Here we have set all physical constants as well as the rest mass of the particles to 
unity, and 

v 

v= - (1.4) 

denotes the relativistic velocity of a particle with momentum v. This system is 
called the relativistic Vlasov-Klein- Gordon system. It is supplemented by initial 
data 

/(0)=/, tt(0)=ui, d t u(0)=u 2 . (1.5) 

The study of this system was initiated in where the existence of global weak 
solutions for initial data satisfying a size restriction was proved. This size restriction 
was necessary because the energy of the system is indefinite so that conservation 
of energy does not lead to a-priory bounds for general data. A natural next step 
in the study of the Vlasov-Klein-Gordon system is the existence theory of classical 
solutions, locally and if possible globally in time. This is the topic of the present 
investigation. 

Another motivation for studying this system of partial differential equations is 
an intrinsically mathematical one. Since the field equation is hyperbolic the system 
resembles the relativistic Vlasov-Maxwell system, for which the quest for global-in- 
time classical general solutions is still open. One might hope that studying related 
systems can help in understanding these open problems more thoroughly. In fact in 
this work we follow the general outline of the existence proof of Glassey and Strauss 
for the Vlasov-Maxwell system. Note, however, that the existence theory of weak 
solutions of the two systems is quite different El E3 • 

The paper proceeds as follows. In Section |2] we prove some a-priori estimates 
necessary for the proof of our main result. These estimates rely on representation 
formulas for the first and second order derivatives of the Klein-Gordon field u, cf. 
Lemmas 12.11 and 12.21 Note that in contrast to the corresponding parts in we 
also need to bound the mixed second order derivatives of the field. In Section[3]wc 
prove our main results, a local-in-time existence and uniqueness result for classical 
solutions and a continuation criterion which says that such solutions can blow up 
in finite time only if the support of / in momentum space becomes unbounded, 
cf. Thms 13.11 and 13.31 In Section 01 we briefly show that the continuation crite- 
rion is indeed satisfied in the one-dimensional situation where so that we 
obtain global classical solutions in that case. Finally some material on the (inho- 
mogeneous) Vlasov equation is collected in an appendix as well as some unpleasant 
technical aspects of the proof of the local existence result, which often have been 
omitted in the treatment of related systems. 

2. A-PRIORI ESTIMATES 

Although our notation is mostly standard or self-explaining we mention the 
following conventions: For a function h = h(t,x,v) or h = h(t,x) we denote for given 
t by h(t) the corresponding function of the remaining variables. For a function 
h depending on the variables x,v we denote its gradient by di x>v )h. By ||.|| p we 
denote the usual L p -norm for pe [l,oo]. The subscript c in function spaces refers 
to compactly supported functions. Sometimes we write z = (x,v) el 3 x K 3 . 
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One main ingredient of our analysis are representation formulas for u and its 
derivatives, which will allow us to establish the necessary a-priori bounds. As point 
of departure we recall that the solution of (|1.2fl is given by 



u(t,x)=uy iom (t,x)+u in h.(t,x), i>0, iei 3 , (2.1) 



where 



u hom (t,x) = J ^ u x {y)dS y --^ J ^ (d x u 1 )(y)-ydS y 
1 f o f , JO 1 f „ , v {Ji(0\'t 



8n 



My)ds y --r [ My) (-M-) \ d y 

\x-y\=t 47r J\x-y\<t V ? / S 



1 f °. f ^ JO 1 f o, ,Jl(0 



— u 2 {y)dS y -— u 2 {y)——dy, 

^tJ\x-y\=t ^J\x-y\<t ? 



is the solution of the homogeneous Klein-Gordon equation with initial data as in 



and £ := y/t 2 - \x-y\ 2 , and 



1 r r < ds , 1 / f , .mo 



Uinh(t,x) = -— / / p(s,y)dS y - h— / / p{s,y)——dyds 



with £ := y/(t — s) 2 — \x — y\ 2 is the solution of the inhomogeneous Klein-Gordon 
equation with vanishing initial data, cf. ^2] or \U\ i Jk denotes the Bessel function. 
To derive formulas for the derivatives of u the differential operators 



S = d t + v-d x , T=-uud t + d x ; lj= X ~ !I 



\x-y\ ' 



which are adapted to our system and have first been introduced in [3] in connection 
with the Vlasov-Maxwell system turn out to be useful. 

Lemma 2.1. (Representation of du) 

Suppose uGC 2 is a solution of the Klein-Gordon equation with p given by 

J for some f GC 1 . Then 



d k u(t, x) = (t, x) + Fl (i, x) + (t, x) + F&t, x) + F* (t, x), k G { 1 , 2, 3, t] 
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where Fq is a linear functional of the initial data only, and 



^J\x- y \<tJ l+u-v ' ' \x-y\' 

F r(t,x) = i- f f a^tu^fit-lx-yly^dvj-^, 



4nJ\x-v\<tJ \x-yf 
f r(.^ x ) = I P(t-\x-y\,y)w k dy, 



\x~y\<t 

1 f* f „r„^MO, 

2 



^Jo J\x-v\<t-s ? 



F&frx) = ~f I —^-jiSmt-lx-yly^dvj^-,.. 



4nJ\x-y\<tJ ' ' \x-y\ 2 ' 



Ffax) =~f I p( Sl y)^{t- S )dyd S . 

47r JO J\x-y\<t-s S 



Here £ = \J(t— s) 2 — \x — y\ 2 and the kernels a k and a* are 

Proof. The proof is a straightforward calculation using 

dx^j^S + Hs^--^)^ dt = -^— {S -v-T) 
.7=1 x 7 

where 5^- = 1 and 5jk = for j ^ fc. Terms involving only the initial data are collected 
in Fq and terms involving Tf are integrated by parts, using the identity (Tf)(t — 
\x-y\,y,v) = d y {f(t-\x-y\,y,v)). □ 



Next we turn to the second order derivatives. 
Lemma 2.2. (Representation of d 2 u) 

Suppose u£C 2 is a solution of the Klein-Gordon equation with p given by 

H1.!A) for some f G C 2 . Then we have for fc,^£{l,2,3,£} 

d H u(t, x) = F ke + Ff s + F& + Ff s + Ff T + F k R % + + F% + F% , 
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where are linear Junctionals of the initial data only, 
Fl e s = -i / / c ke (u,v)(S 2 f)(t-\x-y\,y,v)dv 

4?r J\x-v\<tJ 



x — yY (l+u-v) 2 



1 f [ uktt, , ?.\(Qf\(j. \„ „,i „, „\j., dy \M\ ^ C 



Fts = ^-I jbf{^v){Sm-\x-y\,y,v)dv-^, \bf\< 

47r J\x-v\<tJ \x-y\ z [l+w-v) 6 

F& = ^[ I 'a kl (uj,v)f(t-\x-y\,y,v)dv-^, f a u {u,v)du = Q, 

^J\x-y\<tJ \X-y\ 6 J| w | = i 

Fj& = ±f f d ke (cu,v)(Sf)(t- \x-y\,y,v)dvdy, \d ke \ < 

87r J\x-v\<tJ 



Frt = ^I [e ki (u;,v)f(t-\x- y \,y,v)dv-^- . 

87r J\x-y\<tJ \x-y\ (l+uj-v) 



e kl \< 



1+UJ-V 

C 



2 



± (l _V2zi, 

2(<5fet+'5ft-4t'5ft)-l 



0/7T V2 J|a;-i/|<t 



47T 



y) f-7T^ ( Xfe -Vk)(xt-ye) + {l-fat- S it )^Q)dyds. 

,v-y\<t-s V ? £ 7 

#ere we /icwe set x t —t, y t — s, and u>t = l for notational convenience. 

Proof. The proof is a long and tedious calculation similar to the previous lemma. 
The critical part is, of course, to prove that the kernels a kl appearing in the singular 
TT-terms vanish when integrated over the unit sphere. To give some flavor of the 
respective calculations we outline them in the case of the (most complicated) kernel 
a M for l<k,£<3. We have 

CL — - O " 



(l + \v\ 2 ){l+LU-v) 3 (l + \v\ 2 ) 2 (l+LU- V ) 4 



{l+UJ-v) 2 {l+\v\ 2 )(l+LU-v) 

--: a k£ + a ke + af + af . 



Using d Vt ({ y /l + \v\ 2 + Lu-v)- 2 )=-2(v l + Lu l )(y / l + \v\ 2 + Lu-v)- 3 we find 
f u t duj _ _Vk_Q f duj 



duj 

k i^i(^/r+FPU+w-*)) 3 



V k v £ / — ==T 1777 = -47TW fc W£. 

J\u\=l 



Hence J, i =1 <Xi du = 24TrVkVe. By the identity 



{tJT+W+u-v) 2 (^T+W+^-y) 3j {tJT+W+u-vY 
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we obtain 

[ uif Vk du> 

+-d f (d 1 2ie 

167T 47T 

= ^-vtv k + —6tk, 

which implies Ji^i—j^t du) = — WnVkVi — AnS^i- The remaining two terms can be 
integrated directly to yield 

/ die = —8irv k ve, / diu = AwSki- 

J\io\ = l ~ J\ui\ = l 

Summing up we obtain the desired result 

a H duj = Q. 



The computations for the kernels 

tk _ 2v k (uo-v+\v\ 2 ) 3uj k (uj-v + \v\ 2 ) v k 



(1 + wi) 3 (l + \v\ 2 )(l+u-v) 4 (l + |u| 2 )(l+a;-v) 3 
a tt = - 1 — (3\v\ 4 -{u-v) 2 \v\ 2 -\v\ 2 + 3(u-v) 2 +4lu-v\v\ 2 ) 



(l<fc<3) 



(1+uj-v) 4 

proceed along the same lines. □ 

Note that in contrast to the corresponding analysis of the Vlasov-Maxwell system 
|3j we also need to provide the mixed second order derivatives of the field. The fact 
that the averaging property of corresponding the TT-kernel holds is a pleasant 
surprise in so far as the mixed derivatives do not appear in the field equation. 

We can now use these formulas to derive the necessary estimates on the deriva- 
tives of u. Constants denoted by C are positive and may change their value from 
line to line. If they depend on the initial data this is explicitly mentioned. 

Lemma 2.3. (Estimates on df, du, d 2 u) 

(i) Suppose that f&C 1 is a solution of the Vlasov equation 

Sf{t,x,v) = F(t,x)d v f(t,x,v), f(0,x,v) = f(x,v), fe[0,T[, 
for some F^C 1 . Then we have, with z—{x,v), 

R/(t)||oc<R/||oc+c [ (i+\\d x F OOlWR/OOIMs. 

Jo 

(ii) In addition, assume there exists an increasing function P(t) such that 
f(t,x,v) = for \v\ > P{t), and suppose u is the solution \2.1)) of the Klein- 
Gordon equation for 1 6 [0,T[. Then uGC 2 , and 

|19 ( ,^ )U (t)|| 00 <a(^(l+t) 5 (l+^(t)) 5 +^\i- S )(l + ^( S )) 6 ||^( S )||oe^ , 

where the constant C depends on the norms of the initial data. 
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(iii) On any bounded time interval on which F and P are bounded we have 
||9 ( 2 M)U (i)||oo<cfl+log*( sup H^/OOIIooH 

V 0<s<t Jo J 

where log* (s) = s for < s < 1 and log* (s) = 1 + log(s) for 1 > s, and C de- 
pends on the time interval, the bounds for P and F, and the initial data. 

Proof. The assertion in (i) follows directly from (| A. 3|> below with 5 = and 
G(t,z) = (v,F(t,x)). As to (ii), the estimate of du is a straightforward consequence 
of Lemma 12. II bv using S f = Fd v f and getting rid of the v derivatives via integra- 

o 

tion by parts. Then we estimate each term individually using ||/(f)||oo = ||/ ||oo, cf. 
(IA.2|I with g=0, the estimate 

IT^ 2(1+H2) ' 

and the fact that the ^-integration is only over a finite ball of radius P(t). 

To prove (iii) we use the same procedure with Lemma 12 . 21 replacing Lemma 12. II 
First we need to get rid of the second derivative (i.e., the S 2 f term in F§g). So let 
us first assume Sf € C 1 and consider 

~ l f f f Jtlr. . ti\i a2 r\ n i_ J...AX. d V 



F SS = ^j x J Jc M (co,v)(S 2 f)(t-\x-y\,y,v)dv 



\x-y\ ' 



Observe that 

S 2 f = (SF) -d v f + F-d v (F-d v f) — Fc(v)d x f, c jk (v) = d V] v k . 

Inserting this into Fgg and removing all v derivatives using integration by parts 
we end up with an expression involving only first order derivatives of /, which also 
holds without the SfGC 1 assumption by a standard approximation argument, in 
particular, ueC 2 . 

Now we can estimate each term as before except for the TT-kernel which is the 
most critical one, i.e., 

= ^-l [ a H (u>,v)f(t-\x-y\,y,v)dv-^ 
4tt J\ x -y\ <t J \x-yr 



1 1 



a kl (uj,v)f(s,x+(t — s)u>,v)dvduids, 



where we split the s-integral into two parts, I over [0,t— r] and II over [t — T,t]. 
The first of these can be estimated directly by Clog(f/r). For the second one we 
use the averaging property of a ke to write it as 



11=1 —I hi 



t-T~t~ T J\u\ = \ 



i ke (tu,v) ^f(s,x+ (t— s)u>,v) — f(s 1 x,v)dvdu?j ds. 



Hence by the mean value theorem 

\II\<Ct sup R/OOIIoc- 

t-T<S<t 

Summing up we obtain the estimate 

\F& (*) I < C(\og(t/r) +tN (t)), 0<t<T, 
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where N(t) — sup 0<s<t ||<3 x /(s)||oo- For N(t) < t^ 1 we can choose the optimal value 
T = N(t)~ 1 , otherwise we choose T = t, and this yields 

\F^ T (t)\<C\og*(tN(t)). 

Combining these estimates the remaining claim follows. □ 

3. Existence of classical solutions 

We now have collected all ingredients to show existence and uniqueness of local 
classical solutions of the relativistic Vlasov-Klein-Gordon system. 

Theorem 3.1. (Local existence of classical solutions) 

o 

Let f eC*(R 6 ), Mi £C 3 (R 3 ), v,2 (£C 2 (M. 3 ). Then there exists a unique classical solu- 
tion 

feC\[0,T[xR 6 )), ueC 2 ([0,T[xR 3 ) 

of the relativistic Vlasov-Klein-Gordon system ]1. !)) - ]!. 3\) for someT > 0. satisfying 
the initial conditions \1.5\) . Moreover, 

f(t,x,v) = for \x\ >R+t or \v\ > P{t) 

o o 

where R is determined by f and P is a positive continuous function on [0,T[. 

Proof. We begin with the uniqueness part which relics only on Lemma 12.11 Let 
(f^',u^), (f( 2 \u( 2 >) be two solutions satisfying the same initial conditions, and on 
any compact time interval [0,Tq] on which both solutions exist define / = / W — f 1 - 2 ^ 
and u = u^ — u^ 2 \ Then u satisfies (11.211 with (11.3(1 . and / satisfies 

Sf = d x u^d v f^-d x u^d v f^=d x ud v f^+d x u^d v f. 

Proceeding on a term- by-term basis using the representation from Lemma 12. H and 
replacing the Sf term via the Vlasov equation and integrating by parts we obtain 
the estimate 

||^^(i)||oo <C^* (||/( S )||oo + ||^^(s)|| oc ||/ (1) ( S )|| oc + ||^^< 2 Hs)||oo||./'(^)|| 00 ) 

Using the boundedness of /W and d x u^ this implies 

||3x«(t)||oo<C' / (||/(*)||oo + ||flL«(a)||oo)da. 
Jo 

On the other hand, by (|A.2fl 

f(s,x,v)= f (d x ud v fW)(s,ZW(s,t,x,v))ds, 
Jo 

where Z^ 2 '(s,t,x,v) is the solution of the characteristic system corresponding to 
u( 2 \ Hence we obtain, using the boundedness of d v f^\ 

||/(t)]|oo<C / \\d x u(s)\\oods. 
Jo 

Combining both estimates gives 

\\f(t)\\oo + \\d x u(t)\\oo<C f (||/(s)||oo + ||5»«(a)||oo)da 

Jo 

and Gronwall's lemma implies f(t) = u(t) — 0, proving uniqueness. 
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Next we turn to existence. To this end we set up an iterative scheme and prove its 
convergence to a solution. Let /(°) =/, =u\ and define f^ n > eC 1 ([0,oo[x]R 6 )), 
u {n) gC 2 Qo,oo[xK 3 ) recursively via 

Sf {n) - d x u {n -^d v f [n) = 0, / (n) (0) = / 

and 

d^u {n) -/\u {n) +u {n) = -p {n \ u (Tl) (0)=w°i, d t u {n) (0)=u 2 

where 

p<"> (*,*)= / f^(t,x,v)dv; 



notice that satisfies a support estimate as the one asserted for the solution /, 
but for all £ > and with the increasing function 

P^(i):=sup{M|/( n )(T,a;,t;)^0, 0<r<t, ieR 3 } 

instead of P. 

Step 1 (Uniform bounds on P' n ) and du^): 
By Lemma \'2. 31 (ii) and (|A.2J| we have 

\\d (ttX) u^\t)\\ x <c(il+t)\l+P^) 5 + j\ 
and 

ids, 



p^(t)<p+ ( iia^"- 1 )^)!!, 

Jo 



where P bounds the w-support of the data. 
Define Q( n \t):=max < k < n \\d (t!x) u(V(t)\\ 00 . Then 

g(")(i)<C^(l + t) 5 (l + P^(<)) 5 + ^(i-s)(l + P ( " ) (s)) 6 Q (n) (s)^, 

and by Gronwall's lemma, 

Q( n )(i)<c((l + t) 5 (l + p(")(t)) 5 exp(t 2 (l + p(™)(t)) 6 )). 

If we insert this into the estimate for P^™) we find that 

pM(t)<P+cJ (l + s) 5 (l + p( ,l )(s)) 5 exp(s 2 (l + p( n )( S )) 6 )rfs, t>0, (3.1) 

Hence by induction, P^™' (t) < P(t) , n € N, t £ [0, T[, where P is the maximal solution 
of the integral equation corresponding to l|3.1[l . which exists on some time interval 

o 

[0,T[ whose length is determined by P and the norms of the initial data entering 
the constant C. In addition, d(t, x ) u ^ (t) are bounded in terms of P(t) on that 
interval. 

For the rest of the proof we now argue on a bounded, arbitrary, but fixed time 
interval [0,To] C [0,T[. Constants denoted by C may now depend on To and the 
bounds established in Step 1. 
Step 2 (Uniform bounds on d 2 u^ and d z f^ n '): 
By Lemma \2. 31 (i) and (iii), 

II^^^WIloo < c (n-^* (in- Il^f^^-^C^IIooll^/^C^Iloo)^) , (3.2) 
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\\dl, x) u^(t)\\oo<C (l + log*( sup \\d x f^(t)\\ 0o )+ Wdf^u^is^ds ). 

\ 0<t<T Jo / 

(3.3) 

Applying Gronwall's lemma to (|3.2(l yields 

ll^/^WIIoo^Cexp^C^V + ll^^-^MHoo)^). (3.4) 
Inserting this into 1)3.3(1 we obtain 

II^^WIIoo^C+cAl + ll^^-^^lloo)^, 
JO 

and by induction, 

The bound on ||dJ" ( ")(i)||oo now follows from g2J. 

5iep 5 (Uniform Cauchy property of f (n} (t) and d x u^(t)): 

Introduce 

jm,n _ j(m) _ j (n) ^ u m >" — u ( m ) _ 

and note that 

5/ m '" = (5 x u ( "~ 1) )a„/ m ' n + a a; u m - 1 '™~ 1 5 t ,/ (m) . 
As in the uniqueness part we derive the estimates 



R« m,n (*)||oo<C / (HflxU^^WIIoo + H/^WIloo)^ 

and 

||/ m ' n (t)||oo<C / ||flU« m - 1 ' W - 1 (*)||oods. 

Jo 

Combining these we obtain 

||0 x u m ' n (t)||oo<C / ||5x« m - 1 ' n - 1 (*)||o da, 
Jo 

and hence by induction, 

||&u m '"(0||oo<C^^, fe = min(m,n). 

So d x vS n \t) is a uniform Cauchy sequence, and the same is true for f^ n \t). 
Step 4 (Uniform Cauchy property of d z f^ n \t) and d? t x \U^ n \i)): 
We begin by establishing a bound on \\d? t x s ) u m ' n {t)\\ 00 . Using the representation of 
Lemma |2.2I proceeding on a term- by-term basis and again using a splitting as in 
the uniqueness part we obtain 

II^^^^WIIoo^C^' (||c? ( 2 f , 2;) W r — 1 '"- 1 (3)||oo + ||^r"^(^l|oo) ^- (3.5) 

Here the critical TT-kernel is estimated as in the proof of Lemma 12.31 letting r — > 
t. For the other terms we use the boundedness assertions as well as the Cauchy 
properties already obtained for the lower derivatives. 
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Next we prove that the characteristics converge uniformly. Writing Z^ n \s,t,z) = 
(X( n \s,t,z),V( nS> (s,t,z)), z=(x,v), where we omit the arguments if there is no 
danger of misinterpretation, we find 



— -(x^ n) — x^ m h 

ds 



< |y(")_y(m)| ) 

= |S a; M( ,, - 1 »(xW)-^tt (m " 1) (^ (m) )| 

< |a I !i ( "- 1) (x (n »)-a I « ( "- 1) (i w )| 

+ \d x u^ n - 1 \X {m ^)~d x u {m - 1 \X {m "')\ 

< C\xW-XW\ + 6 m , n . 

Here the expression 5 m ^ n converges to zero if m,n — >oo by the Cauchy property of 
d x u^ and we have used the boundedness of d^.u^ n \ Combining these two esti- 
mates and again using Gronwall's lemma we obtain the claimed convergence of the 
characteristics, i.e., Z^ (s,t,z) converges uniformly for 0<s<T o ; the convergence 
also uniform w.r.t. the parameters t and z. 

Writing Z^ n \s) for Z^ n \s,t,z) the analog of equation (|A.3(1 for iterates and 
vanishing g implies 

dj(t, Z ) = dJ(z^(o))- [ d z fW(s,z^( s )) ^("-^a 1 "^))^, 



and hence 

\d z f m ' n (t,z)\ < \dJ(z^ m \o))-dJ(z^(o))\ 



+ 



d z f (m \s,Z (m) {s)) dlu^-^is^X^Xs)) 

-d z f {n) {s,Z {n \s)) dlu^ n - 1 \s,X { ^(s)) 



ds. 



The first term vanishes in the limit m,n— >oo, and we split the second term into 
the following four parts: 

J' {\d z f^ m \ s ,z^\ s ))-d z f^\ s ,z^\ s ))\ idlu^-v^x^is))] 

+ \d z f {m) {s,Z {n \s))\ |^M( ro - 1 )(s,X( m )(s))-9^( m - 1 )(s,X^(s))| 

+ \d z f im) (s,Z^(s))-d z f^(s,Z^(s))\ Iflgu^-^Cs.JfW^))] 

+ \d z f {n) (s,Z {n \s))\ |a2 w ( m - 1 )( s ,X(")(s))-^ M («-i)( Sj lW( s ))|) ds. (3.6) 



The first two vanish in the limit m,n— >oo by Lemma 13.21 below while the latter 
two may be estimated due to the boundedness of d x u^ and d z f( n \ Hence 

WdJ^it^U^S^ + C [\\\d z f m ' n (s)U + \\d*u m - 1 ' n - 1 (s)\\ 00 )ds, 

Jo 

where again £ mj n— >0 if m,n—>oo. Another application of Gronwall's lemma gives 



||9 z r' m (t,z)|| 0O <<5 m , n + C / \\d x u 



:)2„,m — l.n— 1 



(s) Hoods. 



(3.7) 
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We finally insert the last inequality into equation (|3.5|1 to obtain 

||a? t , a) u m 'ioo<W+c A^)^- 1 '"- 1 ^)!!^ 

Jo 

and, by iteration for each I £ N 

||5f t , x) u™>»||oo<* m>n e c + ^^, 

on [0,To], where B is a bound for ||9? t stt m ||oo- This proves the required Cauchy 
property for ||d 2 t ^"Hoo while the one for ||<9 z /"(s)||oo follows from equation l|3.7jl . 
thereby finishing the proof. □ 

It remains to show that the first two terms in 13.611 indeed vanish if m,n^oo, 
which amounts to showing that d z f^ as well as d? t x ^u^ are uniformly continuous, 
with modulus of continuity uniform in n. One could be tempted to use the mean 
value theorem but suitable estimates on the second order derivatives of /' n ' and 
the third order derivatives of u^ n ' are not available. We begin by defining 

5 n (t,n) := S up[\d z f^(t,z)-d z f ( - n \t,z')\\z,z' eM 6 , \z-z'\<n), 

8 n (t,n) :=sup{|9 ( 2 t ^ )U ("- 1 )(t, a ;)-a ( \ :c)M ("- 1 )(^y)|| a ;,i / eR 3 , |s-y|<»/}. 

Note that both S n and 9 n are bounded uniformly in n. Our desired result now is 

Lemma 3.2. On any time interval [0,Tq] on which the iterates satisfy the bounds 
established in Steps 1 and 2 in the proof of Thm. \3.1\ the following is true: For all 
e > there exist rjo>0 and no G N such that for all n>no 

S n (t,Vo), d n (t,Vo)<£, te[o,T }. 

We defer the rather technical proof of this lemma to Appendix [B] We now 
establish a continuation criterion for the solutions obtained in Thm. 13.11 

Theorem 3.3. (Continuation criterion) Let (f,u) be a solution of the relativistic 
Vlasov- Klein- Gordon system on [0,T[ as in Theorem \S.l\ Then the function 

P(t):=sup{|u| \f(r,x,v)^0, 0<r<i, xeK 3 } 

is bounded on [0,T[ iff ||c? x it(i)|| oc is bounded on [0,T[. Moreover, if T is chosen 
maximally then any of these bounds implies that the solution is global, i.e., T = oo. 

Proof. To prove that a bound on d x u implies a bound on P we integrate the v- 
component of the characteristic system. For the reverse direction we note that we 
can estimate d x u in terms of P exactly as we did for the iterates in Step 1 of the 
proof of Thm. |3~T1 using Lemma \'Z. II 

Assume now that T is chosen maximally, that P is bounded on [0,T[, and T < oo. 
For any to £ [0,T[ we can use the arguments from the proof of Thm. 13.11 to show 
that a solution with data /(to); u(to), dtu(to) prescribed at t — to exists on some 
time interval [^to + ^j except that there is one technical catch here: u(to), dtu(to) 
are not sufficiently regular to qualify as initial data in the context of Thm. 13.11 
But since we already have the solution on [0,io] we can define the iterates used to 
obtain the extended solution as follows: For (f°,u°) we take a global extension of the 
existing solution with the required regularity and with ||/^(t)||oo) I (t) I |oo , 
1 1 9( 4 >a; )U^||oo, ||<9 2 t jjaW^ | |oo bounded in t. Given the (n— l)st iterate we define the 
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nth iterate exactly as before on [0,oo[. Then all these iterates coincide with the 
solution on [0,io], the data term in the formulas for the field, from which the loss 
of derivatives arises, is the one determined by the data at t = 0, and it is straight 
forward to repeat the arguments from the proof of Thm. lXTl to extend the solution 
to some time interval [0,io + <5[. The crucial point now is that the uniform bound 
on the momenta implies that 5 > can be chosen independently of to, cf. and 
the lines that follow. For to close enough to T this contradicts the maximality of 
T. □ 



4. The one-dimensional case 

In this section we illustrate that our continuation criterion from Thm. holds 
and hence classical solutions are global in the one dimensional case, where x,wGM. 
To do so we first need to derive the representation formulas for u and its derivatives 
in this situation. The standard trick to do this is to observe that u{t,x) solves the 
Klein-Gordon equation 1|1.2[) iff w(t,x,£) = u(t,x)e~ 1 ^ solves the wave equation 

$-%-%)w = -e-*p 

with initial data transformed accordingly. This leads to 

u(t 1 x)=u hom (t,x) + u inh (tix), t>0, xeR, 

where Uh olrL (t,x) is the solution of the homogeneous equation and depends only on 
the initial data for u, and 



Uinh(t,x) = -- I / p(s,y)J Q (^ (t-s) 2 -\x-y\ 2 )dyds. 

Jx-(t-s) 



2 

Hence 



1 '* 



d x iii n h(t,x) = —- / (p(s,x+(t—s)) — p(s,x—(t — s)))ds 

* JO 

-Tj / P(s,y) -^7== , %± {x-y)dyda. (4.1) 

IJo Jx-(t-s) ^{t-s) 2 -\x-y\ 2 

Assume now that we have a (local) solution of the Vlasov-Klein-Gordon system in 
the one-dimensional case, with f(t) compactly supported for all t. As before, 

P(t):=sup{|«| |/(r,a;,u)^0, 0<r<t, x£R}. 

Since / is constant along characteristics of the Vlasov equation, 

p(t,x)<2||/|| 00 P(t). 

Integrating the Vlasov equation w.r.t. x and v implies that ||/(t)||i = ||p(i)||i = 

o 

1 1 /||i, and since Ji(£)/£ is bounded on ^>0, 14.1|l implies that 

||5^(i)IU<C7(l+t 2 + i J P(t)). 
Integrating the w-component of the characteristic system implies that 

rt 

P(t)<P+C(l + t) 3 + C sP(s)ds (4.2) 

so that by Gronwall's lemma P can not blow up on any bounded time interval. 

Indeed, these estimates could be repeated for iterates defined as in the proof of 
Thm. 18.11 Controlling first order derivatives of / ( -™- ) and second order ones of u^ n ' 
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would be much easier than in the three dimensional should be obvious from 

comparing Eqn. 14.1JI with Lemma \'2. II Hence: 

o 

Theorem 4.1. Let f €C'(1 2 ), u\ eCf (M), u 2 £C$(M). Then there exists a unique 
classical solution 

/eC^^oofxl 2 ), ueC 2 ([0,oo[xR) 

of the one- dimensional relativistic Vlasov- Klein- Gordon system, satisfying the ini- 
tial conditions j 1.5)) . If P denotes the solution of the integral equation corresponding 
to £P$ then 

f(t,x,v) = for \x\>R+t or \v\>P{t) 

o o 

with R determined by f. 

Appendix A. Some facts on the Vlasov equation 

In this appendix we collect some useful facts on the (inhomogeneous) Vlasov 
equation for easy reference. As before, we combine x and v to one variable z = (x,v), 
and we consider the initial value problem 

d t f(t,z) + G(t,z)d z f(t,z)=g(t,z), f(0,z) = f(z), 

o 

where / , G g £ C 1 and G is such that the solutions of the corresponding characteristic 
system 

z(s) = G(s,z(s)) 

exist on the time interval on which G is defined. Denote by si— > Z(s,t,z) the solution 
corresponding to the initial condition Z(t,t,z) = z and recall that 

Z(s,t,z) = Z(s,r,Z(r,t,z)) 

and 

d t Z(s,t,z)+d z Z(s,t,z)G(t,z) = 0; 

the second equation follows from the first by differentiating z = Z(s,t,Z(t,s,ZQ)) 
with respect to t and then choosing zq = Z{s,t,z). The partial derivative d z Z{s,t,z) 
satisfies the first variational equation 

d z Z(s,t,z) = {d z G){s,Z{s,t,z))d z Z(s > t > z), 

with d z Z(s,s,z) =1 the unit matrix, or equivalently, 

d z Z(s,t,z)=I+ J (d z G)(r,Z(r,t,z))d z Z(r,t,z)dr. 

In addition, there is also the less obvious equation 

d z Z(s,t,z)=I+J d z Z(s,r,Z(r,s,z))(d z G)(r,Z(r,t,z))dr, (A.l) 

which holds because the right hand side solves the first variational equation. How- 
ever, note that the integrands in the last two equations are not equal. 

We now apply these results to the Vlasov equation. Clearly the solution is given 

by 

f(t,z)=f(Z(Q,t,z))+ f g(s,Z(s,t,z))ds. (A.2) 
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Moreover, d z f(t,z) exists and satisfies 

d z f(t,z) = (dJ)(Z(0,t,z))- f (d z f){s,Z{s,t,z)(d z G)(s,Z(s,t,z))ds 



(d z g)(s,Z(s,t,z)ds. (A.3) 

This can be seen by a straightforward calculation using equation l|A.l(l . We should 
remark that the usual way of deriving this equation is by differentiating the Vlasov 
equation with respect to z and viewing the resulting equation as an inhomogeneous 

o 

Vlasov equation for d z f. This requires f,g&C which was not necessary here. 
Appendix B. Proof of Lemma EOl 

We split the proof into several steps, following the approach used in ^3 pp. 
78-90] in the case of the Vlasov-Maxwell system. Throughout we argue on a time 
interval [0,T ] on which the iterates satisfy all the bounds established in Steps 1 
and 2 of the proof of Thm. 13.11 

Lemma B.l. (Estimate on 9 n in terms of 9 n -\, 5 n -\) 

There exists C > such that Ve > 3% > Vr? € [0,770] V< € [0,T ] Vn e N, 

9 n (t,ri)<£ + C (v + J (6 n -i(s,ri) + 5 n -i(s,ri))ds\ . 

Proof. Using once more the representation formulas of Lemma 12.21 we proceed as 
in Steps 1 and 2 of the proof of Theorem 13. II □ 

Lemma B.2. (Estimate on 8 n in terms of 9 n ) 

There exists C > such that Ve > 3t]q > V77 G [0,770] V< € [0,T ] Vn e N, 

6 n (t,rj)<E+C [ 9 n (s,r))ds. 



To prove Lemma lB.2l we need some additional information on the derivatives of 
the characteristics which is provided by the following lemma. 

Lemma B.3. (Estimates on the derivatives of the characteristics) 

(i) There exists C>0 such thatVs,tE[0,T } Vz,z'gR 6 VnGN 

\Z^(s,t,z)-Z^(s,t,z')\<C\z-z'\ 

(ii) There exists C>0 such thatVs,te[0,T }, s<t, Vz^'eM 6 VneN 
\d z Z^(s,t,z)-d z Z^(s,t,z')\ 

<C\z-z'\+J^ l^u'"" 1 ) (r,!'"' (r,t, z)) - S^'^ 1 ) (t, jt"' (r,t, z')) | ^. 

Proof. Part (i) follows immediately from equation IA. lit . As to (ii), we only treat 
the terms involving ^-derivatives; the c^-terms may be estimated in the same fash- 
ion. So we start with the term \d Xi X^(s,t,z) — d Xi X^ (s,t,z')\, l<i<3. By an 
elementary calculation, cf. PP- 83-85], we obtain 

\d s d Xi X^(s,t,z)-d s d Xi X^(s,t,z')\ 
<2\d Xi V^(s,t,z)-d Xi V^(s,t,z')\ 

+ 5\d Xi V^(s,t,z)-d Xi V^(s,t,z')\ \V^{s,t,z)-V^{s,t,z')\. (B.l) 



16 



M. KUNZINGER, G. REIN, R. STEINBAUER, AND G. TESCHL 



Using 

d s d Xi V^(s,t,z) = -d Xi d x u^- 1 \s,X^(s,t,z)) d Xi X^(s,t,z) 

we obtain 

\d s d Xi V^(s,t,z)-d s d Xi V^(s,t,z')\ 

^dlu^^s^Xs^-dlu^is^Xs^z')^ \d Xi X^\s,t,z)\ 
+ \d 2 x u^- 1 Hs,X^(s,t,z'))\ \d Xi X^\s,t,z)-d Xi X^\s,t,z')\. (B.2) 
Using (i) and equations IB. 1(1 and ljB.21) we obtain 
\d Xi Z( n \s > t > z)-d Xi Z^(s,t,z')\ 

<c(\z-z'\ + J \d Xl Z^(T,t,z)-d Xt Z^(T,t,z r )\dr 

+ ^|9> (n " 1) (T,X(")(r,<,z))-^ u (n-i) (T;X (™) (Tjij/)) | (ir y 

Another appeal to Gronwall's lemma completes the proof of the lemma. □ 

Proof of Lemma \B.2\ We only treat the <9„/l(n)-terms; the c^/^-terms can be 
estimated analogously. Again writing z— (x,v) respectively z' = (y,w) we find 

\d v fW(t,z)-d v fW(t,z')\ 

< \dj(z<ri(p,t,z))\ \d v z^(o,t,z)-d v z^(o,t,z')\ 

+ \d v Z^(0,t,z')\ \dJ(ZW(0,t,z))-dJ(ZW(0,t,z'))\. (B.3) 

The second term in the above estimate is bounded by e for large n, t£ [0,T ] and 
\z — z'\ suitably small by the fact that d z Z^ is uniformly (in n) bounded, the 
uniform continuity of the data d z f and Lemma [B. 31 (i). 

By Lemma IB. 31 (ii) and the fact that the characteristics are uniformly bounded 
the first term in equation (|B.3|I may be estimated by 

Cx\z-z'\+Cxf dlu^^ (t,X^ (r,t,z))- d 2 x u^^ {t,X^ (t, t, z')) 



dr. 

Denote by C 2 the maximum of C\ and C from Lemma fB. 31 (i). Set C = (\2C2~} + 
l)Ci where [r] denotes the smallest integer bigger or equal to r. We claim that 
this constant verifies the assertion. Indeed, let e>0, and 776 [0,770] with 770:= 
e/(2C\ + l). Let z,z'eR 6 with \z-z'\ <r], nGN and r,tG [0,T ]. Then by Lemma 
R3l(i) 

\Z ( - n XT,t,z)~Z ( - n \T,t,z')\<2C 2 7 1 . 
Together with the fact that 9 n (T,krj) = k9 n (T,rj) for all fceN this gives 

|^(-*- 1 )(r ) lW(r ) M))-3 t 2 U («- 1 )(r 1 XW(r,t )2 '))| 

< 6 n ( T ,2C 2V ) <9 n (T,(l+\2C 2 -\) V ) < (1+ {2C 2 ])9 n (T,ri). 

Summing up we obtain 

|a„/< n > (t, z) - 0„/ W (t,z')\<e + C f 6 n (r, 77) dr. 

Jo 

□ 
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Proof of Lemma Vi.'A Combining Lemma IB. II and Lemma IB. 21 we obtain: 



Ve>0 3r ?0 >OVr ? G[0,%], te[0,T ], ne N: 6 n (t,rj) < C + J q O n -i{s,V)ds 



which by iteration shows that 



Ve>0 3?7o>0 3n eN Vn>n , te[0,T] : 6 n (t,r] )<e. 



The claim concerning S n now immediately follows from Lemma IB. 21 □ 
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